Exact results for the Kondo screening cloud of two helical liquids 
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We analyze the screening of a magnetic quantum dot with spin i coupled to two helical liquids. 
Interestingly, we find two qualitatively different sets of Toulouse points, i.e., nontrivial parameters 
for which we can solve the two channel Kondo model exactly. This enables us to calculate the 
temperature and voltage dependent Kondo screening cloud, which develops oscillations for an applied 
spin voltage /i s . Such a spin voltage can be conveniently applied by a charge bias in a four-terminal 
helical liquid setup. 



o 



-a 
c 

o 
o 



00 

o 
[> 

l> 

O 

(N 



X 



PACS numbers: 72.10.Fk, 73.43.Nq, 75.76.+J 

Introduction. The crucial ingredient of Kondo physics 
is the coupling of a localized spin degree of freedom, 
often represented by a spin on a quantum dot (QD), to 
a spin bath. In recent years, peculiarities relating to a 
plethora of realizations of the spin bath have been in- 
vestigated 0^3- One of the archetypal phenomena in 
Kondo physics which is still subject of active research is 
the characteristic screening of the QD-spin: while theo- 
rists predict a macroscopically extended screening which 
has been coined the Kondo cloud, experimental confir- 
mation of this unique correlation is still lacking 0, 
Problems related to the direct detection of the Kondo 
cloud are the high frequency at which the spin of the QD 
is flipping and the principal inability to directly mea- 
sure correlations between the QD-spin and the lead-spin 
without decisively perturbing the tunneling region. He- 
lical li quids, w hich have recently been theoretically pre- 
dicted [10l4l2| and experimentally discovere d 11311 a t the 
edge of the quantum spin Hall insulator [l(J llll , fea- 
ture two modes of excitations with opposite spins moving 
into opposite directions. The quantized conductance of 
G = e 2 /h of a single helical edge state is topologically 
protected against backscattering by time-reversal sym- 
metry. A magnetic quantum dot coupled to helical liq- 
uids is one of the simplest nontrivial perturbations which 
allows elastic backscattering assisted by a spin flip of the 
QD 0]. Along these lines, former publications @, 
have mainly focused on the effect of a magnetic impu- 
rity on the conductance of the helical liquid and not on 
the screening of the localized spin. However, attaching 
helical liquids to a magnetic QD also offers unique oppor- 
tunities to investigate spin-dependent scattering off the 
QD: by carrying away spin resolved information about 
the Kondo QD in distinct directions, correlations become 
measurable away from the QD avoiding the experimental 
necessity of locally perturbing the tunneling region. 

We investigate an S — \ QD that is coupled to two he- 
lical liquids in a two channel Kondo model [16]. The rele- 
vance of this model is corroborated by a renormalization 
group study in Ref. [17]. In this Letter, we provide ex- 
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FIG. 1: (Color online) Two helical liquids of a quantum spin 
Hall insulator are coupled via an anti-dot. Properly gated, 
the dot resembles a magnetic impurity with S = = . Each side 
is brought out of equilibrium by an applied voltage. 



act results for the local screening and the Kondo screen- 
ing cloud. To this end, we solve the model by restrict- 
ing the parameters to configurations for which an exten- 
sion of the method of Emery and Kivelson [18j allows us 
to map the full interacting Hamiltonian to a quadratic 
one. We call these configurations Toulouse points fol- 
lowing the terminology of the noninteracting case. In- 
terestingly, we find two qualitatively distinct sets A and 
B of Toulouse points; set A contains the Toulouse point 
of the two channel Kondo model for noninteracting leads 
[lif , while set B resembles a one channel Toulouse point 
at both Luttinger parameters g = 



15| . Determin- 



ing the local screening and the spatially extended Kondo 
screening cloud, we are able to demonstrate a different 
phenomenology for the novel interacting set of Toulouse 
points. For case A, the spin is always fully screened, 
however, case B, i.e. the interacting regime, provides 
the full range from no screening to exact screening. The 
Kondo cloud, defined as the correlations between QD- 
spin and lead-spin, obeys an asymptotic decay for large 
x that is quadratic at zero temperature. At finite temper- 
ature T — l/(/cg/3), the decay becomes exponential after 
a length scale £t ~ hv/3, which was also seen by Borda 
[81 for the one channel Kondo model. In contrast to the 



noninteracting one channel case, we find a In 2 (a;) instead 
of an x _1 divergence for small x. Furthermore, we ob- 
serve that an applied spin flavor voltage fi x = fi t + fi hl 
whereby "flavor" denotes top (t) or bottom(b) (seeFigJT}, 
and a spin voltage jj, s = fit — fib have distinct effects on 
the Kondo cloud. While the former acts as an artificial 
magnetic field which decreases the extent of the Kondo 
cloud significantly, the latter induces spatial oscillations 
of the Kondo cloud. Hence, the Kondo cloud can be eas- 
ily manipulated in our four-terminal setup. 

Model. A schematic setup is drawn in Fig[T] where 
the magnetic QD is realized by a properly gated anti- 
dot. Employing an anti-dot rather than a adatom in 
connection with helical liquids has the advantage that 
the anti-dot conserves the same pseudo spin, described 
in Refs.flil]. as the helical liquids. This is important 
as the normal spin is not a good quantum number for 
experimentally realized helical liquids due to spin orbit 
coupling. 

The modeling Hamiltonian is the two channel Kondo 
model for helical liquids, given by the sum of the three 
terms 

H a = dxl ^2v ayF ¥(x) a:a (aid x )^(x) a ^ 
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for each side o G {t,b}. Here, the operators ^> a ,a and r 
represent the fermions of the leads and the spin of the 
QD respectively, v a ^F is the Fermi velocity, the interac- 
tion strengths within the leads are denoted by (72/4,0 and 
the coupling constants J° with J% = J% = determine 
the strength of the coupling of the leads to the QD. Fol- 
lowing Ref. 20], the helical liquids are bosonized using 



/27ra 



the bosonization identity ^ a>cr (x) 
and diagonalized by introducing the fields 
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with the Luttinger liquid parameters g a and u a , which 
depend on v a ,F and 52/4,0- This leads to 

Ua = J dx {dxVaA x )) 2 1 

-J^F^ei^^r- +H.c, 
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TABLE I: Toulouse points of the two channel Kondo model 
for helical liquids. There are two disconnected sets A and B 
of Toulouse points. The parametrization uses v G (0, 00) and 
qe (0,tt/2). 



We furthermore include the chemical potentials fi a by the 
nonequilibrium operator 15l. l21| 



y= Aa(^.,t-^« r 



a£{t,b} 



(4) 



where M a ,a is the total number operator of the fermion 
^> a ,a- In principal, independent chemical potentials at 
each terminal depicted in FigQ] could be considered. 
However, as the Kondo interaction is a spin interac- 
tion, the two potential configurations that are only able 
to alter charge densities are not relevant for the Kondo 
physics. Therefore, we do not consider them here. 

Method. The aim is to find the Toulouse points 
of the Hamiltonian, i.e. points in the parameter 
space for which there is a mapping of the Hamiltonian 
from Eq.© to a quadratic one following the ideas of 
Refs.[3, |22|, [23|]: (i) applying an Emery- Kivelson ro- 
tation U — exp(i ^ a X ai+ (p ai+ (0)T z ), (ii) transforming 
ip a ,± orthonormally to fields fa = Y, a ,a M o,(a,<j) Va,a, (in) 
rcfcrmionizing. This technique has also been recently em- 
ployed to analyze the Kondo problem for a single helical 
liquid in Ref.Q. 

In our setup, we obtain all Toulouse points by restrict- 
ing J* to the value such that H.KZ.a cancels the Emery- 
Kivelson rotation of H a and secondly imposing that all 
vertex operators after step (ii) take a refermionizable 
form. Here we distinguish two possibilities. Within case 
A both vertex operators take the form e ±l ^^°\ while 
for case B the form of one vertex operator deviates to 
e ±i0 2 (o\ where 2 and 04 are linearly independent. A 
representation of all Toulouse points is given in Tab|U 
Case A is characterized by gt + gb = 2 and contains the 
noninteracting case. In contrast, the novel case B obeys 
9t + <?b = 1 and intrinsically relies on interactions. In 
particular, it contains gt = .% = \ for which refermion- 
ization is known to be a promising method in similar 
models 15, 2(| and possesses two solvable values of J^ h 
for each solvable configuration of gt/b- 

The grand canonical operators of the resulting reso- 
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nant level models are 18 1 
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with the local pseudofermion c = G\t~ , for case A, and P ear f «r g t = g b = \ for which J z = J z = J z 



Hb = -(m* - Ms)A/2 - {Hx + VsWi + H 



J^l(0)c + J^l(0)c + H.c 



(6) 



with c 



r for case -B. The fermions are determined by 
-j W — -/27r a Gje^i ( x > , where we exploited the treat- 
ment of Klein factors Gj and number operators Afj in 
Refs.[20l. [23| for the cases A and _B, respectively. In 
both cases, the noninter acting Hamiltonian is given by 
Ho = v^2j J dx ^>j(x)(—id x ) i Sj(x). The main difference 
between the Hamiltonians is the absence of double cre- 
ation terms in Hb- In particular, Hb shows no signs of 
Majorana fermions, which are central in the treatment 
of Ha 18j] • This is a direct consequence of the intro- 
duction of two linearly independent fermionic fields in 
case B. Interestingly, for both cases, the dependence on 
the Luttinger parameters on both sides is excluded from 
the Hamiltonian and incorporated in the transformation 
between the original fermions ^> a and As will be 

seen below, this gives rise to simple algebraic relations 
between observables at different Toulouse points. 

The Hamiltonians of Eqs.© and (JSJ) are solvable via 
a variety of techniques, see, e.g. Refs.|2ll |23{ . Our re- 



sults are based on infinite order perturbation theory in 
the Keldysh formalism. Conveniently, we exploit a treat- 
ment introduced in Ref.[l5| and replace the voltage /i x 
by an artificial magnetic field H alt at the QD for the 
Hamiltonian of Eq. ([5]) 



-/x x A/" 2 -> frc'c =: ff art /is 5Qd(t 2 



(7) 



with H art — ^i x /{^b 5qd)- This substitution leaves 
all thermodynamic expectation values considered in this 
paper invariant and has an appealingly simple physical 
foundation given in Ref. [15| . Notably, for a Lande factor 
5qd on the order of 10, a spin flavor voltage of IV gen- 
erates an artificial magnetic field on the order of 10 3 T. 

Results. We focus now on the z-screening of the spin 
r of the QD. The contributions to this quantity are the 
local screening at x = 0, determined by the locally bound 
spin in the leads, and the spatially extended Kondo cloud. 
The local screening Sq is given by 
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E 




with p spin .a = |(p ,t - Pa,i), e = 0+, and the "=F"-sign 
corresponds to the two possible branches of Toulouse 
points in case B described in Tab|U Hence, for all 
Toulouse points of case A the spin is locally exactly 
screened, case B however provides Toulouse points rang- 
ing from zero to exact screening. Both extremes ap- 

At 

this Toulouse point, the local screening is more meaning- 
fully described in the form Sq = — J z /(4ttv)(t z ). Con- 
sequently, the local screening vanishes for the Toulouse 
point with J z = and reaches exact screening for 
J z = An v. 

The Kondo screening cloud Xa( x ) ^ s represented by the 
spatially resolved correlations between the QD-spin and 
the lead-spin 25] 

xl{ x ,9t,b) = (<5p S pin,aO)£T 2 ), (9) 

with 8 A = A — (A). We find that in our case, all Kondo 
clouds are derivable from two generic Kondo clouds. One 
stems from the noninteracting Toulouse point in case A 
and the other from the gt = gb = h Toulouse point in case 
B. We denote them by Xai x , 1) an d Xa( x > \) respectively. 
For case A, the general Kondo cloud is related to the 
noninteracting one by 

X z a {x,9t,b) =gaX z a {x,l) (10) 
with a € {t,b}. Furthermore, the sum rule 

-Xt{x,9t,b) + — Xb(<B,0t,b) =0. (11) 
9t 9b 

connects the Kondo clouds of both sides. In combination, 
Eqs. fTTl) and (fTU)) imply that there is no Kondo cloud in 
case A if j£- = Jg } i.e. in the case of symmetric coupling. 
Interestingly, Eq. (ITT1) also shows that the Kondo correla- 
tions of different sides have different signs. Thereby, the 
side of positive correlations is determined by the weaker 
in plane coupling | | . 

Similarly, for case B, the Kondo clouds for different 
Luttinger parameters are related by 

Xb( x ,9t,b) = V9b ((VSb - V9t)xl( x , |) 

+( v / ^ + VsDxbM))' ( 12 ) 

and the Kondo cloud on the t-side follows by exchanging 
b <-> t in the above formula. A direct consequence of 
Eq. (fT2|) is the independence of the total Kondo cloud on 
the Luttinger parameters of both sides 

X z t(x,9t,b)+X^(x,gt,b) = Xt(x, \) + X z ,(x, !)■ (13) 

In equilibrium, the analytic expressions for the generic 
Kondo clouds read 



X t z (x,l) =sgn(J, 
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FIG. 2: (Color online) The decay of the generic Kondo 
clouds in equilibrium for different inverse temperatures f3* — 
kgT^/gP on a double-logarithmic scale. For finite temper- 
ature, the Kondo cloud decays exponentially after a length 
scale £t ~ hvfi. The relevant limits of the cloud are depicted 
as dashed curves, 7 is the Euler-Mascheroni constant, (a) 
Noninteracting Kondo cloud. The temperature of reference 
is T% ■ Both Kondo temperatures T^ 2 determine crossover 
length scales of the cloud at £[y 2 = to/(isTy 2 ). (b) In- 
teracting Kondo cloud rescaled by T t x . The temperature of 
reference is Tg . The curve is universal for all and only 
the length scale determined by Tg is relevant. 



where $(x,P,T K ) = $(e"^t, 1, \ + ^j3k B T K ) with the 
Hurwitz-Lerch transcendent 3>(z, s, a) — Y^ 00 n 7—% — w, 



-1/2 



are the Kondo temperatures of our 

"t^ = i67rfciija • equilibrium Kondo 
cloud is shown in Fig(3] Concerning case A, given in 



system [21| and T t 



Fig |2(a)[ we chose Tf = y / TfJf to be the temper- 
ature of reference. Both Kondo length scales = 

hv/(kBTy 2 ) appear in the shape of the Kondo cloud 
by inducing a crossover to different asymptotic behav- 
iors. The leading divergence for small x is ln 2 (x) and the 
asymptotic behavior for large a; is a quadratic decay. For 
finite temperature the Kondo cloud decays exponentially 
in x after a length scale £t ~ fvvff similar to Ref.Q. 

For case B, the temperature of reference is conve- 
niently chosen to be T§ = Tf + T 2 K '. Rescaled by T t K , 



FIG. 3: (Color online) The generic Kondo clouds rescaled by 
x 2 at zero temperature in nonequilibrium. (a) Noninteract- 
ing Kondo cloud. An applied spin voltage fi s — kgTjf/g/j,* 
induces exponentially decaying oscillations of the frequency 
fi s /(hv) in the Kondo cloud and locally changes its sign, (b) 
Interacting Kondo cloud. Additional permanent oscillations 
of the frequency 2/j, 3 /(hv) appear. 



the generic Kondo cloud of side t is a universal curve for 
all coupling parameters. In contrast to case A, the only 
relevant length scale is determined by £g = hv / (fcsTg ). 
The asymptotic behavior of the depicted curves in Figj2] 
coincides in the large x limit and differs only by the con- 
stant ln 2 (r i Rr /T 2 ic )/(47r) 2 for small x. For Tf = Tf, 
they become equivalent. 

Applied voltages fi x and /z s alter the Kondo cloud dis- 
tinctly. The spin flavor voltage \i x acts as a magnetic field 
and strongly shrinks the Kondo cloud in both classes of 
Toulouse points. This reflects the sensitivity of Kondo 
physics to unhindered spin flips of the QD. An increas- 
ing fj, s instead induces oscillations in the Kondo cloud de- 
picted in Fig(3J For case A, an oscillation of the frequency 
Us/ (hv) changes the sign of the Kondo cloud locally. This 
oscillation decays exponentially and the Kondo cloud be- 
comes monotonic again for large x. For case J5, the uni- 
versality of the Kondo cloud at a specific Tg is destroyed. 
In addition to the known oscillations of case A, case B 
shows oscillations of the frequency 2fi s /(Hv), which decay 
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quadratically in x. Nevertheless, the Kondo cloud stays 
negative for all x. The oscillations can be interpreted 
as the impact of Friedel oscillations on the Kondo cloud. 
While Friedel oscillations in helical liquids are suppressed 
in the presence of ordinary scatterers, the spin-flip scat- 
tering off the magnetic QD to the oppositely moving 
channel generates an interference of the wave functions 
in both channels. 

Conclusions. We couple a magnetic QD with spin \ to 
two helical liquids at the edges of a quantum spin Hall in- 
sulator. For the two channel Kondo model, we calculate 
the Toulouse points at which the Hamiltonian is mapped 
to a quadratic one. In contrast to non-interacting leads, 
the Toulouse points cover all possible coupling parame- 
ters to the QD at the cost of restricting the interaction 
strengths and Fermi velocity of the helical liquids. The 
Toulouse points separate into two qualitatively different 
sets. The first is related to the known Toulouse points 
of the non-interacting case. The second appears exclu- 
sively in the case of interactions. We determine the local 
screening and the Kondo screening cloud and calculate 
the signature of applied voltages. 
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